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SUMMARY 


PROBLEM 

Determine  time  evolving  spectrum  statistics  for  acoustic  signals  that  have  rever- 
berated from  the  ocean  floor. 

RESULTS 

A theoretical  investigation  based  upon  the  ray-wave  synthesis  of  the  time  evolv- 
ing spectrum  of  ocean  floor  reverberation  has  been  completed.  This  study  is  applicable 
to  the  case  of  moving  platforms,  arbitrary  beam  patterns  and  a varying  sound  speed 
structure. 

RECOMMENDATION 

Apply  the  results  of  this  study  to  practical  ocean  acoustics  information  transfer 
problems. 
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1.  INTRODUCTION 


Acoustic  echo  ranging  signals  and  acoustic  communication  signals  that  are  received 
in  the  ocean  frequently  are  distorted  by  reverberation  generated  by  the  interaction  of  the 
transmitted  signal  with  the  ocean  bottom.  Additionally,  if  the  transmitter  and  receiver 
platforms  are  moving,  then  the  Doppler  shifts  in  frequency,  induced  by  the  platform  move- 
ments, further  complicate  the  reverberation  effects  on  the  received  signal.  This  gives  added 
cause  to  consider  reverberation  to  be  an  undesirable  noise  component.  In  the  usual  case  of  an 
assumed  known  environment  the  reverberation  statistics  can  be  determined  to  enable  the 
transmitter  and  receiver  to  be  designed  to  combat  this  noise. 

The  problem  considered  here  is  an  arbitrary  but  narrow  band  signal  that  is  trans- 
mitted by  one  moving  platform  and  received  by  another  moving  platform.  The  platforms 
are  assumed  to  have  arbitrary  transmit  and  receive  beam  patterns,  respectively.  The  sound 
speed  structure  is  assumed  to  be  spatially  varying  in  a simple  manner  and  the  direct  path 
signal  generally  follows  a bending  ray  path  to  the  receiver.  In  addition  to  the  direct  path 
signal,  the  receiver  intercepts  the  reverberation  generated  by  the  interaction  between  the 
transmitted  signal  and  the  ocean  bottom.  It  is  assumed  that  the  random  bottom  has  a 
gaussian  distribution.  The  time  evolving  spectnim  of  the  field  is  determined,  which  re- 
veals interesting  frequency-spatial  relations  and  provides  useful  estimates  of  the  acoustic 
field. 

This  development  is  restricted  to  the  case  of  bottom  scattering  to  avoid  the  mathe- 
matical complexity  imposed  by  a moving  surface.  However,  it  is  possible  that  the  bottom 
reverberation  results  can  be  applied  to  the  case  of  surface  reverberation,  with  certain  quali- 
fications concerning  Doppler  spreading  of  frequencies  caused  by  surface  movement.  The 
effects  of  volume  scattering  are  ignored  also  as  this  reverberation  exhibits  more  compli- 
cated frequency-spatial  relations  than  boundary  scattering.  This  is  not  a serious  omission 
when  bottom  reverberation  is  present,  because  volume  reverberation  usually  is  much  weak- 
er than  bottom  reverberation. 

The  analysis  leading  to  these  results  appears  in  Chapter  2 where  the  developed 
basic  field  expression  represents  the  propagation  and  scattering  mechanisms  described 
earlier.  This  development  is  based  upon  a generalization  of  the  Kirchoff-Helmholtz  inte- 
gral equation  ^ >2  in  which  the  Green’s  function  kernels  have  been  replaced  by  ray  theory 
kernels.  The  surface  field  in  this  integral  is  approximated  by  a ray-Bom  approximation, 
i.e.,  the  surface  field  is  equated  to  the  direct  path  ray  field.  The  surface  field  derivative 
is  approximated  by  assuming  local  specular  reflection  with  shading,  which  is  a slight  gen- 
eralization of  the  Kirchoff  approximation.  This  development  is  patterned  after  that  of 
Eckart^,  who  used  the  conventional  Helmholtz  integral  equation  and  approximated  the 
boundary  conditions  via  the  Bom  and  Kirchoff  approximations.  Eckart’s  boundary  assump- 
tions were  later  extended  by  Horton  and  Muir^  and  Horton,  Mitchell  and  Barnard^  to 
include  a shading  function.  Additionally,  his  method  of  approximating  the  integral  was 
later  extended  by  Gulin. h 

In  each  of  these  earlier  developments,  as  well  as  in  more  recent  ones,  the  primary 
interest  has  been  centered  on  the  case  of  specular  reflection  and  its  attendant  reverbera- 
tion, i.e.,  the  scattered  return  from  the  first  few  Freznel  zones.  An  exposition  of  much 
of  this  earlier  work  is  given  in  Beckmann  and  Spizzichino^  and  a more  recent  survey  is 
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given  by  Fortuin.8  The  specular  problem  is  discussed  here,  as  well  as  the  reverberation 
that  follows  the  specular  component,  i.e.,  from  the  higher  order  Freznel  zones.  While 
this  component  of  reverberation  has  not  been  thoroughly  researched  one  of  the  early 
investigations  of  this  reverberation  was  conducted  by  Halley.9>KM  1 

In  Chapter  3 the  mean  and  covariance  function  of  both  components  of  rever- 
beration are  found.  These  results  give  the  statistical  properties  of  the  time  dependent 
field  that  are  needed  for  the  determination  in  Chapter  4 of  the  time  evolving  “gener- 
alized instantaneous  power  spectrum”  (GIPS)^  and  a]so  of  the  spectral  variance.  The 
GIPS  of  the  reverberation  from  the  higher  order  Freznel  zones  reveals  interesting  fre- 
quency-spatial relations  and  represents  a generalization  of  Halley’s  earlier  work. 

2.  SCATTERING  BY  THE  OCEAN  BOTTOM 
THE  TRANSMITTED  SIGNAL 

A transmitter  and  receiver  are  assumed  to  be  traveling,  respectively,  at  arbitrary 
but  constant  velocities  vj  and  vj^,  as  indicated  in  figure  1 . The  transmitter  emits  a 
signal 


s(t)  = m(t)  exp(-i  cogt) 

(la) 

m(t)  = \ f M(a>)  exp(-i  cot)  dco 

Z7T./-00 

(lb) 

that  propagates  to  the  receiver  along  the  direct  ray  path  and  also  along  ray  paths  that 
scatter  from  the  bottom.  The  emitted  signal  is  Doppler  shifted  by  the  motion  of  the 
transmitter  and  amplitude  weighted  by  the  transmitter  beam  pattern  and  the  direct  path 
signal,  as  well  as  the  signal  that  scatters  from  the  bottom,  retains  this  information.  The 
received  signal  is  Doppler  shifted  further  by  the  motion  of  the  receiver  and  weighted  by 
the  receiver  beam.  An  expression  is  developed  for  the  received  signal  in  which  the  effect 
of  each  of  these  mechanisms  is  apparent.  For  mathematical  convenience  the  case  of  scat- 
tering by  the  moving  surface  is  ignored;  i.e.,  only  the  case  of  scattering  by  a stationary 
bottom  is  considered.  However,  the  bottom  is  assumed  to  be  acoustically  penetrable  with 
scattering  coefficient  R(fp)  and  it  is  assumed  that  for  the  appropriate  scattering  coefficient 
the  development  approximately  applies  to  the  case  of  surface  scattering. 

THE  HARMONIC  FIELD 

The  field  in  the  medium  due  to  a harmonic  transmission  at  frequency  co  is  given 

byl>2 

<t>u>(r)  = *5^(  /dv  [gco(£)  + t(|,  r)] 

V 

+ / ds  n • [0w(?)V-V  <Ml)l  } X(l,r)  (2a) 
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Figure  1.  Ray  geometry  from  transmitter  to  receiver. 


(2e) 


7(£>r)  = V2  A(|,r)/A(|,r) 

where  denotes  a vector  dot  product,  the  time  dependence  in  equation  2a,  exp(-i  cot),  is 
omitted  for  convenience  and  the  ray  kernel  is  defined  in  equations  2b,  c and  d.  The  quantity 
n(=  cq/c ) is  the  index  of  refraction,  where  c is  the  sound  speed  and  cq  is  an  arbitrary 
reference  sound  speed,  the  quantity  k(=  co/cq)  is  the  wave  number;  the  integrals  in  equation 
2c  and  equation  2b  are  over  ray  paths  between  the  points  r and  f. 

For  the  case  of  a constant  sound  speed  structure,  the  phase  expression  in  equation  2c 
becomes  the  range,  r,  the  amplitude  expression  in  equation  2b  becomes  the  reciprocal  of  range, 
!/r,  and  the  function  y vanishes,  which  implies  that  there  is  no  volume  scattering.  It  is 
assumed  that  the  sound  speed  is  varying  sufficiently  slowly  so  the  volume  scattering  term 
in  equation  2a  can  be  ignored.  This  formalism  is  strictly  valid  only  for  the  case  of  a nonmov- 
ing transmitter  and  receiver  but  for  velocities  v that  are  small  compared  to  the  sound 
speed,  the  associated  error  terms  are  of  the  order  Ivl/c  and  can  be  ignored. 

Initially  the  transmission  is  considered  from  one  element  of  the  transmitter  ar- 
ray to  one  element  of  the  receiver  array.  The  transmitter  source  term  is  given  by 

gg(£)  = 47r6($-fj)  Fp^d-p)  (3) 

where  5(*)  denotes  the  delta  function  and  Fp^dj)  denotes  the  frequency  dependent  com- 
plex transmitter  transfer  function.  The  field  at  an  element  of  the  receiver  array  due  to 
this  transmission  is 


^co^R’  M = 1Tw(^T)  x(*T’ 

+ Ttt  /ds  " * &r>  V - V p. 

X(fp>  Ir)  (4) 


where  the  transmitter,  receiver  and  surface  scattering  points  are  denoted  explicitly  by  £7. 

and  £p.  This  is  an  integral  equation  which  can  be  evaluated  to  first  order  by  a ray- 
Bom  approximation  and  a modified  Kirchoff  approximation. 


Boundary  Conditions 

The  field  on  the  scattering  surface  can  be  expressed  as  the  sum  of  the  incident 
and  reflected  fields 

0a>dpdj)  = ^codpdj)  + dp>£r)  * (5) 

The  incident  field  is  approximated  by  the  shadowed  direct  path  signal  while  the  reflected 
field  is  approximated  by  the  associated  specular  reflection. 

The  amplitude  and  phase  of  the  specular  component  will  be  the  same  as  the  inci- 
dent component  except  as  they  are  modified  by  the  reflection  process  and  the  total  field 
on  the  surface  can  be  written 

*w(Wt>  = n(ip,|x)  rTcj(fT)  x(lT,IP)(i  + R(IP))  (6) 

where  IT  is  the  shadowing  function  which  is  unity  on  the  ensonified  areas  and  zero  on  the 
shadowed  areas. 

Similarly,  the  field  gradient  term  is  given  by 

n ’ ^ 0co(£p’It^  = ” ' V[0*(£p>fj)  + 0r(fp,£j)l 

= i k n(£p,£x)  r'XtiJ(£'p)  x(£X;£p) 

X n • [V^'dpix)  + R(fp)  1 (?) 

where  the  amplitude  gradient  terms  are  much  less  than  the  phase  gradient  terms  and  are 
ignored.  Since  the  phase  expression  in  equation  2c  is  symmetrical,  the  gradient  of  the  phase 
of  the  incident  component  is  approximated  from  the  direct  path  phase 

V^'dp.t-p)  = VWMp) 

= v^dpiT) 

= nC^pl'edpJy) 

where 't\|p,|x)  is  the  unit  vector  that  is  directed  along  the  ray  from  the  transmitter  to 
the  surface  at  fp. 

The  gradient  of  the  phase  of  the  specular  component  immediately  above  the  scat- 
tering plane  can  be  found  from  the  condition  that  its  vector  wave  number  has  the  same 
projection  in  the  scattering  surface  as  the  incident  component 

n x ( V i//' + V i//r)  = 0 • (9) 

Since  at  the  scattering  surface  the  vector  wave  numbers  have  the  same  magnitude 

\W\=WV\  (10) 
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their  projections  perpendicular  to  the  scattering  surface  are  also  the  same 

n • (V<//*  +V^0  = 0.  (11) 

With  this  result  the  field  gradient  expression  becomes 

n • v |j)  = i k xi(|p,|j)  rXcj(IT>  X(IT.Ip) 

Xn  •'SXIP,IT)  (1-R«P))  (12) 

where  the  index  of  refraction  is  ignored  since  it  is  of  the  order  of  unity. 

The  gradient  of  the  ray  kernel  is  similarly  given  by 

n V x(fP,fT)  = * ^ X(IpJr)  n ■'^(Ip>^r)  (13) 

where 'Stlp.lp)  is  the  unit  vector  that  is  directed  along  the  ray  from  the  receiver  to  the 
surface  at  fp. 

Collecting  terms  the  expression  for  the  field  at  the  receiver  element  becomes 

0co(Ir>It)  = x(It’£r) 

+ ^ /ds  Hdp.lp)  rTwdx)  xdjdp)  X(^pJr) 

X n • [/SdP,|R)(  1 +R(|P) ) - ^|P,IT)(  1 -R(|p) ) 1 • 

(14) 

By  Huygen’s  principle  (see  appendix  A)  the  surface  integral  over  just  the  direct  path  term 
vanishes,  i.e.,  the  integral  vanishes  for  R = 0,  and  this  expression  is  simplified  to  the  form 

^codR-M  = fTcjdx)  X(Ix-Ir) 

+ f ds  TT(|p,|j)  rTcjdT)  xdxdp)  X(Ip'Ir)  R(Ip) ) 

X n • feX£p,£R)  + 'eX£pJ£x)  1 (15) 

which  is  a first  order  approximation  of  the  field  at  the  receiver  elements  at  £p  due  to  a 
harmonic  transmission  from  the  transmitter  element  at  fx- 

Beaniformed  Signals 

The  beamformed  receiver  output  caused  by  the  beamformed  harmonic  signal  trans- 
mission can  be  obtained  by  summing  over  the  transmitter  and  receiver  apertures.  For  this 
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calculation  it  is  assumed  that  the  ranges  between  the  arrays  and  the  bottom  are  large  com- 
pared to  both  the  array  sizes  and  bottom  elevations  and  approximations  are  obtained  for 
the  phases  and  amplitudes. 

The  phase  of  the  direct  path  signal  is  expressed 


r m 

\K£t>£r)  = du  • ne(uJR) 


(16) 


where  u denotes  an  arbitrary  path  between  £R  and  |j.  This  path  is  chosen  via  the  ray 
path  to  the  acoustic  center  of  the  transmitter, 

| 

«Mr)  = T du  • ne(u,|R)  + «//(?T,|R)  . (17) 

rT 

Since  the  integrand  is  nearly  constant  for  ranges  that  are  large  compared  to  the  array  sizes, 
the  phase  approximated  is  • 


\K£t>£r)  - eTR  ' (£7  " rT^  ^rT’^R^  (18) 

where  eyR(=  e(fj,rR)  s e(ry,£R))  is  the  unit  transmitter  polar  vector  directed  along  the 
ray  from  the  receiver  acoustic  center  to  the  transmitter  acoustic  center  and  n = 1.  Identi- 
cal arugments  yield  the  expression 

</'(Ir,£t)  = ®RT  ' (£r“^r)  + ’/'(^R^T^  (19) 


where  eRj  is  the  unit  receiver  polar  vector  directed  along  the  ray  from  the  transmitter 
acoustic  center  to  the  receiver  acoustic  center.  By  symmetry  of  the  phase  function 
equation  1 8 becomes 


'/'(It-Sr)  ~ &pR  • (fc'p-rj)  + eRj  • (£R-rR)  + ^(r-p,rR)  • (20) 

The  amplitude  of  the  direct  path  signal  is  approximated 

A($'P’£r)  = Mrj,rR)  (21) 

since  the  field  expression  is  relatively  insensitive  to  amplitude  variations. 

Similarly,  for  the  phase  and  amplitude  expressions  within  the  integral 
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Mfrrip)  “ %>  ‘ (^T"rT^  +^PT  * k ^rP^  + ’KrT,rp)  (22a) 

A(£y,£p)  = A(r-p,rp)  (22b) 

<Mfp ,fR)  = ^PR  • k r)(r?)  +%p  • (fR-rR)  + i//( rp,rR)  (22c) 

A(|p,|R)  = A(rp,rR)  (22d) 


where  r\  denotes  the  surface  elevation,  ?p  denotes  the  projection  of  the  surface  point  £p 
onto  the  plane  of  the  mean  surface  and  the  vectors'©^  and^py  are  the  unit  polar  vectors 
at  ru  and  ?p,  respectively,  directed  along  the  rays  from  ?p  and  ?u. 

Collecting  the  linearized  phase  expression  and  substituting  into  equation  15,  along 
with  the  variable  change 

ds  = dxdy/n  • Jc  (23a) 

n = (z-T?)/ 1 V(z-r?)  | , (23b) 

multiplying  the  complex  receiver  transfer  function  rRa)(|R)  and  integrating  over  the 
apertures 


^co(rR’rT) 

“ P^-Str)  Poo(-£rt)  A(rj,rR)  exp(i  k i|/(rj,rR)) 

+-^  1 dxdy  II(rp,rT)  A(rT,rp)  A(rp,rR)  R(rp) 

X P^I-eyp)  Paj(“®Rp)  V(z-r?(rp)  )(epj+epR) 
X exp(i  k[i//(rT,rp)  + ^(rp,rR) 

+ 77(rp)  k • (eV-  ,/s  ' 1 ) 

(24) 

Pg/^Tu) 

= fdij  fxojdj)  exP(i  k e-pu  • (|y-rj) ) 

(25a) 

Pco^Ru* 

= /<%  rRo;dR)  exp(i  k eRu  ■ (IR-rR) ) 

(25b) 

where  the  functions  II  and  R were  redefined  to  be  functions  of  the  mean  plane  coordinate 
and  the  function  II  is  approximated  as  a function  of  the  transmitter  acoustic  center.  The 
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functions  p^C-ejy)  and  P^l^Ry)  represent  the  complex  beam  patterns  of  the  trans- 
mitter and  receiver  beams,  respectively,  and  the  beam  “look”  direction  is  defined  to  be 
opposed  to  the  unit  vector  arguments,  i.e.,  back  along  the  rays.  For  the  arguments  ejR 
and  eRy  the  complex  beam  weighting  is  given  in  the  direction  of  the  ray  that  propagates 
between  the  transmitter  and  receiver.  The  beam  patterns  in  the  surface  integral  with 
arguments  eyp  and  e^p  represent  the  cross  sections  of  the  complex  beam  patterns  of 
the  transmitter  and  receiver  that  lie  in  the  scattering  surface. 

THE  PULSED  FIELD 

Given  the  above  expression  for  the  beamformed  output  due  to  a beamformed  har- 
monic transmission,  the  analogous  expression  for  pulses  is  obtained  by  a Fourier  trans- 
formation, provided  that  the  spatial  length  of  the  pulse  is  long  compared  to  the  array 
dimensions  and  bottom  elevations.  The  Fourier  transform  of  the  complex  envelope  of 
the  transmitted  signal  is  denoted  in  equation  lb  by  M(tu)  and  the  Fourier  transform  of  the 
analytic  signal  is  Mfco-cog).  where  cjq  is  the  angular  carrier  frequency.  The  received  sig- 
nal is  given  by 


In  this  Fourier  transformation  the  beam  patterns  are  assumed  to  be  constant  over  the  band 
of  the  signal,  the  integral  coefficient  is  approximated  by  kg(=  cjq/cq),  which  is  valid  for 
narrow  band  signals,  and  the  small  order  term  in  the  complex  envelope  is  ignored.  This 
expression  represents  the  field  in  the  medium  in  terms  of  a broad  class  of  signals,  beam 
patterns,  sound  speed  structures  and  bottoms  for  the  “frozen”  case  and  in  the  following 
section  the  effects  of  platform  movement  are  included. 


0(rT,rR.O  = P(-®TR.)  P^RT^  A^T’^r) 


X m(t-\[/(rT,rR)/co) 

X exp[-i  oJo(t-’//(rT^R)/cO)l 


X p(-e  pp)  p(-eRp)  V (z-rj(rp) ) • (epY+^PR) 
X m(t-('/'(rY>rp)  + \Krp,rR)  )/cq) 

X exp[-i  cuo(t-(^(rY,rp)  + \Krp,rR)  )/cq) 


+ i kQ  7?(rp)  k • (epx+epR)  1 • 


(26) 
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PLATFORM  MOVEMENT 


For  the  case  of  a moving  transmitter  and  receiver  the  platform  position  vectors, 
rj(t)  and  rR(t),  are  assumed  to  be  functions  of  time  while  the  motion  of  the  medium  is 
ignored  for  simplicity.  The  surface  motion  has  the  effect  of  Doppler  broadening  of 
spectral  components  while  the  internal  motions  of  the  mediums  have  the  effect  of  Dop- 
pler broadening  and  biasing  of  spectral  components.  Both  effects  are  assumed  to  be 
negligible  compared  to  the  Doppler  shifts  due  to  the  transmitter  and  receiver  motions. 

For  the  moving  transmitter  and  receiver,  the  direct  path  phase  is  given  by 


\Krj(  r),rR(t) ) 


rT(r) 

R(r) 


du  • n e(u,rR(t)) 


(27a) 


r = t - ^(rj(r),rR(t))/co  (27b) 

where  r is  the  time  at  which  the  received  signal  was  transmitted.  If  the  travel  distance 
of  the  platforms  is  short  compared  to  the  range  between  platforms  during  the  interval 
of  a pulse  it  can  be  linearized 

^(rT(r),rR(t))  = f T du  • n e(u,rR(t))  + vKrj(O)  ,rR(t) ) 

rT(0) 


= e(rT(0)  ,;R(t) ) • vTr  + iJ/(rT( 0)  ,?R(t) ) • (28a) 

The  second  term  is  found  from  the  symmetrical  relation 


- - f ?R(t)  - 

\KrR(t)  ,rT(0) ) = Jh  du  • n e(u,rT(0))  + i//(rR(tD),rT(0) 

rR(tD) 


= e^R^o)  >rT(°^  ' ^(t-to)  + ^R^D^  ) 

(28b) 

tD  = i//(rT(0)  ,rR(tD)  )/cq  (28c) 

where  tjj  is  the  arrival  time  via  the  direct  path  of  the  transmission  at  r=0.  Terms  col- 
lected are 
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(29a) 


> >rR(t)  ) - CjR  • VjT  + e^-p  • v R(t— tj-j) 


+ \KrT(0),rR(tD)) 


T 


(wd) 


i-CRT  • vR/c0 
l+eTR  • vT/cQ 


(29b) 


where  the  unit  polar  vectors  opR  (=  e(rj(0)  »rR(tp) ) s e(r-p(0),rR(t)))  and 
eRT(=  ®(rR(tp>)  »fj(0) ) are  generalized  to  denote  the  ray  directions  at  the  transmitter 
position  when  the  transmission  starts  and  on  the  receiver  position  when  the  reception  begins, 
respectively.  This  occurs  when  the  receiver  is  at  its  position  of  initial  reception  or  the 
transmitter  at  its  position  of  initial  transmission. 

Similarly,  for  the  phase  in  the  scattering  integral 


\KrT(r),rp)  + i//(rp,rR(t) ) 


r rT(r)  - --  r rP 

= / du  • n e(u,rp)  + du  • n e(u,rR,(t) ) 

J r„(t) 


rP  *R 

r = t - [ ^(r-p( r ) ,rp)  + i//(rp,rR(t) ) ] Icq 


(30a) 


(30b) 


where  r is  the  time  at  which  the  signal  received  via  rp  was  transmitted.  These  expressions 
can  be  approximated  by 


\!/(rj(r  ) ,rp)  + i^(rp,rR(t)) 

®TP  ' vTt  + ®RP  * vR(t-tPR)  + ^(rT^®)  >rp)  + *//(rP’rR(tPR) ) 

Ola) 


r = 


(t-tPR) 


1-gRP  ‘ vr/c0 
1+ejp  • vT/co 


(31b) 


tpR  = l\Kr*p(0)  ,rp)  + ^(rp,rR(tPR) ) 1 /cq  • (31c) 

where  the  subscript  R on  the  time  delay,  tpR,  emphasizes  that  the  receiver  position  at  the 
time  of  reception  is  different  for  different  scattering  points.  A first  order  expansion  about 
the  receiver  position  at  the  onset  of  the  bottom  return,  rR(tp),  i.c.,  from  the  specular 
point,  gives 
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Old) 


\Krp.rR(tpR))  - ^(rp,rR(tp))  + $Rp  • VR(tpR-tp) 

where  the  unit  vector  eRp(=  e(r(tp),rp) ) is  defined  to  be  the  unit  vector  at  the  receiver 
that  is  directed  along  the  ray  from  rp  to  the  receiver  location  at  the  onset  of  the  bottom 
return.  Substituting  in  equation  31  gives 

i//(rT(r ) ,rp)  + \[/(rp,rR(t) ) 

= eTP  • vjr  + feRp  • vR(t-tp) 

+ i/zCr-pfO)  ,rp)  + </'(rp>rR(tp) ) 

1_6RP  ' vR/c0 

r = (t-tp)  — — — 

1+eTp  ■ vT/c0 

tp  = [ i^(rT(0),  fp)  + i//(rp,rR(tp)  ) ] / cQ 

where  the  two  receiver  displacement  terms  combine  to  give  a total  displacement. 

The  amplitude  expressions  are  approximated 

A(rj(r)  .fR(t) ) = A(rj(0)  ,rR(tp) ) (33a) 

A(rT(r ) ,rp)  A(rp,rR(t) ) = A(rj(0)  ,rp)  A(rp,rR(tp) ) (33b) 

since  the  direct  path  term  and  scattering  integral  are  relatively  insensitive  to  small  amplitude 
variations. 

Collecting  results  equation  26  becomes 

0(rR»rT »0  = P(~£tr)  P(-eRj)  A(rj,rR) 

x m(t-tD)  exp[-i  w0  a D^t-tD^ 
ik0  f 

+ I dxdy  ll(rp,rj)  A(r-p,rp)  A(rp,rR)  R(rp) 
x p(-ejp)  p(-eRp)  V(z-rj(rp) ) • (epy+epR) 
x m(t-tp)  exp[-i  coq  Op(t-tp)  + ikg  flOp)*^!  (34a) 


(32a) 

(32b) 

(32c) 
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where 


aD  - l 1 - (ejR  • vT  + ej^j  • vr)/cq] 
ap  = [ 1 - (8Tp  • vT  + eRp  • vr)/cq] 

vz  ~ * ' (gPT  + ^PR) 


(34b) 


(34c) 


(34d) 


and  the  small  order  terms  in  the  complex  envelope  are  ignored.  These  equations  very  nearly 
represent  the  final  field  equations  but  they  can  be  simplified  further 

THE  RECEIVED  SIGNAL 

In  the  preceding  calculations  it  was  convenient  to  consider  the  unit  ray  vectors  to 
be  directed  away  from  their  corresponding  source  points.  In  the  following  it  will  be  con- 
venient to  reverse  their  direction,  i.e.,  to  negate  all  unit  ray  vectors.  Then,  whereas  the 
vector  euv  was  defined  to  be  the  vector  at  ru  directed  away  from  rv  along  their  adjoining 
ray,  i.e.,  along  the  ray  that  appropriately  leads  or  follows  the  platform  movement,  it  is 
now  defined  to  be  the  vector  at  ru  directed  toward  rv  along  their  adjoining  ray,  i.e.,  along 
the  ray  that  appropriately  leads  or  follows  the  platform  movement. 

The  amplitude  expressions  were  retained  in  their  general  functional  form  through- 
out the  preceding  calculations  but  for  the  following  they  will  be  approximated  by  the 
spherical  spreading  law 


where  r is  the  transmitter-receiver  range  and  rj  and  rR  are  the  distances  between  the  trans- 
mitter and  receiver  and  the  mean  scattering  plane,  respectively. 

With  the  variable  change  t - t+tp  equation  34  becomes 


A(rT»rR)  - 1/r 

A(r-p,rp)  A(rp,rR)  = l/r-prR 


(35a) 


(35b) 


<t>( t)  = yP(ejR)  P(eRj)  m(t)  exp(-i  tug  «d 


X V (z-7j(rp) ) • (cpx+®PR^ 

X m(t-rp)  exp[-i  ojq  ap(t-r  p)  -ikQ  t?( ?p)^z]  (36a) 


- 1 1 + (ejR  • vj  + eRj  • vR)/c0] 


(36b) 
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Ofp  = [1  + (eTR  • vT  + eRp  • vR)/c0]  (36c) 

vz  = k • (%  + epR)  (36d) 

t p = tp  - tj)  (36e) 

P(®Tu)  = J d^T  ^Tco^T^  exP(-'  k0  ®Tu  ' ) (360 

P(^Ru)  ~ J d^R  rRcj^R)  exP("‘  k0  ®Ru  ' (£r_ir)  ) (368) 

which  are  our  principle  field  equations. 


The  functions  II,  R and  r?  are  assumed  to  be  independent  random  variables.  Their 
distributions  are  unknown  but  by  the  central  limit  theorem,  the  integral  term  in  equation 
36a  has  approximately  a Gaussian  distribution.  Consequently,  the  total  field  is  a finite 
duration  narrow  band  signal,  i.e.,  the  direct  path  signal,  in  a nonstationary  Gaussian  ran- 
dom signal.  The  field  is  completely  defined  statistically  by  its  mean  value  and  correlation 
function  and  the  Fourier  transform  of  its  correlation  function  gives  the  G1PS  of  the  pro- 
cess. These  functions  are  calculated  in  Chapters  3 and  4. 


3.  THE  ACOUSTIC  FIELD  STATISTICS 


THE  MEAN  FIELD 

The  mean  value  of  the  field  is  the  sum  of  the  direct  path  term  and  the  mean  of  the 
bottom  scattering  term 

<0>  = 0o+<^l>  <37a> 

<t>o  = ' P(ejR)  P(crt)  m(t)  exP(-i  ^0  (37b 

<*<>  = 

X p(eTp)  p(eRp)  Vz  m(t-r  p) 

X exp[-i  cjq  otp(t-rp)  ] Q(k0  vz)  (37c) 

where  the  mean  value  of  the  bottom  height  is  zero  and  the  characteristic  function  is  de- 
fined 
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Q(k0  vz)  = <exp(-i  kQ  17  t-z)>  • 

(38) 

If  the  bottom  elevations  satisfy  a Gaussian  distribution 

p(i?)  = (27rh2)_1/2  exp(-r?~/2h-) 

(39a) 

> 

(39b) 

then 

Q(k0  vz)  = exp[y(k0  h vz)2) 

(40) 

where  h is  the  mis  bottom  height. 

At  any  time  t the  contributions  to  the  reflected  signal  come 
portions  of  the  bottom  specified  by 

from  the  ensonified 

0 < (t-r  p)  < T 

(41a) 

T P = *P  ‘ lD 

(41b) 

tp  = [ iMrT(0)  ,rp)  + i//(rp,rR(tp) ) ] /cQ 

(41c) 

= )/cq  • 

(4  Id) 

The  term  tp  represents  the  propagation  time  to  the  intersection  of  an  ellipsoid-like  surface, 
with  the  transmitter  and  receiver  at  its  focii,  and  the  mean  scattering  plane,  while  the  term 
tp  represents  the  direct  path  propagation  time  between  the  transmitter  and  receiver. 

The  onset  of  the  direct  path  signal  occurs  at  the  time 

t0  - 0 • 

(42a) 

Denoting  the  specular  point  in  the  mean  plane  r§p  the  condition  in 
that  the  onset  of  the  specular  scattering  occurs  at  the  time 

equation  41a  shows 

II 

CO 

(42b) 

and  terminates  at  the  time 

t2  =rSP  + T • 

(42c) 
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During  the  time  interval  tj  < t < t2  the  reflected  signal  comes  from  a region  within  an 
ellipse-like  boundary  centered  about  the  specular  point.  For  the  times  t > t2  the  re- 
ceived signal  comes  from  a region  confined  between  two  ellipse-like  boundaries  centered 
about  the  specular  point.  Denoting  the  vertical  projections  of  the  transmitter  and  re- 
ceiver positions  onto  the  mean  plane  by  rjp  and  r^p,  the  outer  ellipse-like  boundary 
passes  under  the  transmitter  and  receiver  at  the  respective  timesr-pp  andrj^p,  while 
the  inner  ellipse-like  boundary  passes  under  the  transmitter  and  receiver  at  the  respec- 
time  times  (rjp  + T)  and  (r^p  + T).  The  average  times  of  passage  are  defined  as 

tj  = (t  jp  + 7"  Rp)/2  (42d) 

*2  = ('TP+'RP>/2  + T ’ (42e) 

The  ordering  of  these  events  is  dependent  upon  the  pulse  length,  but  for  signals  that  are 
short  compared  to  the  range  between  the  platforms,  say,  the  events  are  ordered 

*0  ^ ^ *2  ^ ^ *2  * (420 

During  the  time  interval  tj  < t < t2  the  reflected  signal  is  denoted  specular  scattering  and 
it  comes  exclusively  from  the  region  between  the  transmitter  and  receiver,  i.e.,  the  “in- 
terior” region.  For  wave  lengths  that  are  short  compared  to  the  ensonified  region,  the 
exponential  function  in  equation  37c  oscillates  through  a number  of  “lowest  order 
Freznel  zones”  as7p  varies  over  the  interior  region.  During  the  time  interval  t^  < t < °° 
the  reflected  signal  is  denoted  reverberant  scattering  and  it  comes  exclusively  from  the 
complimentary  “exterior”  region  which  contains  the  higher  order  Freznel  zones. 

Over  the  interior  region,  the  radial  vectors  ejp  and  ej^p  are  approximately  con- 
stant, which  implies  that  the  range  of  Doppler  shifts  is  small  and  the  integrand  is  approx- 
imately constant.  Over  the  exterior  region,  the  radial  vectors  range  over  their  maximum 
domain  of  2ir  steradians  and,  for  the  case  of  moving  platforms  and  nonconstant  beam 
patterns,  the  range  of  Doppler  shifts  is  maximized  and  the  integrand  exhibits  its  maxi- 
mum variations.  In  the  following  it  is  convenient  to  consider  just  the  cases  of  specular 
and  reverberant  scattering  for  the  short  pulse  case  and  infer  more  complicated  cases  from 
these  results. 

Specular  Scattering 

Specular  scattering  is  largely  from  a region  around  the  specular  point  rgp  and  the 
specular  scattering  integral  is  approximated  in  equation  37c 
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<0S>  = 


ko<nspXfesp> 

4?ri  rTSP  rRSp 


pCejsp)  P^RSP^  m(t-TSp) 


X f dxpdyp  uz  exp[-  y (vz  kg  h)2] 

X exp[-i  cog  “SP^  ‘ rP-*  1 (43) 

where  the  functions  extracted  from  the  integral  are  slowly  varying  in  the  interior  region 
and  are  approximated  by  their  values  at  the  specular  point.  The  remaining  integral  is 
approximated  by  integrating  a first  order  expansion  of  the  exponential  function  over  the 
lowest  order  Freznel  zones. 

The  time  delay  between  the  transmitter  and  receiver  via  the  scattering  point  rp 
is 


tp  — [ vKr-p(O),  rp)  + \Krp,rR(tp)  )]  / Cq  (44) 

«• 

where  rR(tp)  denotes  the  receiver  position  at  the  onset  of  reception  of  the  scattered  sig- 
nal. This  expression  can  be  expanded  about  the  specular  point  rgp  by 


*P 


n e(u,rR(tp))  du 


lc0 


+ [ ^(ryfO)  ,r$p)  + '/,(rsP’rR(tp) ) ) /cq  (45) 


where  ep-p  and  epR  are  the  unit  vectors  at  the  scattering  point  that  are  directed  along  the 
rays  to  the  transmitter  position  when  the  transmission  starts  and  the  receiver  position  when 
the  reception  begins,  respectively. 

First  order  expansions  of  the  vector  products  can  be  developed  from  figure  2. 

At  the  specular  point 


(eSPT  + eSPR)  * i 


-cosO  + cosO 


= 0 


(46a) 
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Figure  2.  Ray  vectors  near  specular  point. 
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CO  > 


(46b) 


(e$pj  + egpR)  -f  = cosrr/2  + cos7t/2 
= 0 

(e<jpT  + egpR)  • k = cos(ir/2-0)  + cos(7r/2-0) 
= 2 sinO 


(46c) 


where  0 is  the  grazing  angle.  Perturbing  the  reflection  point  to  rp,  which  is  assumed  to 
be  a small  displacement  compared  to  the  range,  then  the  first  order  expansions 


(%p  +^pr)  •'T'  = - cos(0  + eT  ) + cos(0  + eR  ) 


:Rx' 


- sin0  (eTx  - eRx) 


= - (xp-xSp)  sin20(2/R) 

(SpT  +/^pr)  *t  = cos(jt/2  + eT  ) + cos(7r/2  + eR  ) 


(46d) 


Ty-' 

= - (eTy  + eRy) 


where 


= - (yP  - ySP)  (2/R) 

(©)yp  +^pr)  ' k = cos(7T/2-0  + eTz)  + cos(tt/2-0  + eRz) 
= 2 sin0 

R = 2 rTSprRSp/  (rTSP  + rRSP)  • 


(46e) 


(460 


(46g) 


The  bounds  on  the  incremental  angles  are  denoted  by  double  subscripts,  e.g.,  ejx  is  the 
incremental  angle  between  the  ray  to  the  transmitter  and  the  x-axis 

Substituting  these  expansions  about  the  specular  points  into  equation  43,  the  inte- 
gral factor  becomes 

1 = ^zSP  ^Pl-^O'zSP  k0  h)2l  exPH  w0  “SP^t-rSP^  ) 

X f dxpdyp  exp[i  -Q^SP  (sin20(xp-xsp)2  + (yp-ysp)2)  ] • 

(47) 
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With  the  variable  changes 


ll  = 


^O^SP  sin0(xp-Xgp) ; t2  = ^yP  " ySP*  (48a,b) 


this  integral  can  be  cast  jnto  the  form  of  a product  of  Freznel  integrals 


1 = ‘'zSP  exPt’T^zSP  k0  h)“l  exPl_i  w0  “SP(t*TSP)  1 
X(kOaSPsin0)  /dtl  eXp0  *!>  / dt2  exp(i  <2* 


ilZzSP  R — r J_r„  l k\2 


kO«SPs^0 


exp  I-  — (^zsp  k0  h)Zl  exP  ^0  aSP^t-tSP^  1 (49) 


Substituting  back  into  the  specular  scattering  integral,  then 


<0 SR>  ~ <nSP><RSP>  exPt-y(vzSP  k0  h)"l 


X p(ey§p)  P(crsP^  ” TSP^ 

X exp[-i  coq  <*sp(t-rSp)  ^ ^rTSP  + rRSP^ 
[ 1 + (Cygp  • nj  + Cj^gp  ■ rij^)  / Cq] 

2 sin0  • 


aSP 


^zSP 


(50a) 

(50b) 

(50c) 


This  is  a principal  result  showing  that  the  mean  reflection  from  the  near  surface  is  a 
specular  reflected  signal  attenuated  by  shading  absorption  and  surface  roughness.  The 
starting  and  stopping  times  of  this  signal  are  given  in  equation  42b, c. 

Reverberant  Scattering 

For  the  case  of  reverberant  scattering  it  is  convenient  to  represent  the  signal’s 
complex  envelope  by  the  integral 


m(t-rp)  = I dt ' m(t'  ) pe(t-rp-t') 


(51a) 
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0<t'  < e 


Pe(f)  = 1/e 

= 0 • elsewhere  (51b) 

Substituting  into  the  expression  for  the  mean  value  of  the  reflected  field  in  equation  37c, 


<*  R>  = 


Lim  ^0 
e-*0  47ri 


dpyd0y 

rTrS 


X <lI(rp,rT)><R(rp)  > v exp[-^(t>z  kQ  h)2] 

X pf(t-rp-t ' ) exp[-i  cuq  Qfp(t-rp)  ] (52) 

where  pj,  are  the  cylindrical  coordinates  in  the  mean  scattering  plane  about  the  trans- 
mitter. 

The  only  contribution  to  the  inner  integral  at  a time  t satisfying  equation  52 
comes  from  an  ellipse-like  band  defined  by  the  condition 

0 < (t-rp-t')  < e • (53) 

The  radial  coordinate  integration  over  this  band  can  be  easily  performed.  For  a fixed  time 
t,  angular  coordinate  and  parameter  t'  the  upper  and  lower  radial  coordinate  limits  are 
of  the  form 


pT,u  " 

py(t,  0y,  t ' ) 

(54a) 

^T,L  = 

Py(t,  0ys  t 1 + e) 

(54b) 

Expanding  the  lower  limit  in  a Taylor  series  about  the  parameter  t' 


4.  ta; 

PJX  ~ pT,u  + I ^T 


and  the  integral  in  equation  52  becomes 


f(t>0'p,t  ) 


(54c) 
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<^R>  " 


f 2ir 

dt'  m(t')  / 

J° 


d0nr  ^n(rp,ijV>  <R(rp)> 

TjTj 


x p(%-p)  p(%p)  Vz  exp l-j(vz  ko  h)2] 


X exp(-i  co0  apt'  ) (-dpT/dt'  ) ](t  t,  } (55) 

where  the  integrand  is  evaluated  at  the  point  on  the  ellipse-like  band  (t,  Oj,  t').  Since 
the  quantities  in  brackets  are  relatively  insensitive  to  variations  in  t'  and  become  inde- 
pendent of  t'  for  t •*  00,  the  integrals  can  be  approximately  uncoupled 


<^R>  " 4^  / . d0T 


2ir 

0 


[^L<»<Vt)X^p» 


X p(%p)  p(%P)  vz  exp[- ^(^  kQ  h)z] 


X (-d  Pr/dt1  ) / dt1  m(t ')  exD(-icjnapt) 

(t,  0T)  /o 

^ (56) 

Since  the  complex  envelope  is  bandlimited  below  the  frequency  wgap 


0R  = 0 (57) 

which  shows  that  the  signal  scattered  from  the  exterior  region  is  incoherent  and  has  no 
specular  component. 

THE  CORRELATION  FUNCTION 

The  correlation  function  for  the  total  field  is  given  by 

<0(t+r)0  * (t)>  = 0Q(t+r)  0Q*(t)  + 0Q(t+r)  < <f>  j *(t)> 

+ 0O(t)  <<Mt  + T>>  + <0i(t+T)0*(t)> 


(58) 
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where  the  mean  was  computed  in  the  mean  field  section.  The  correlation  function  of  the 
reflected  field  is  given  from  equation  36a  by 


C(t+r,t)  = <<A](t+r)  0|*(t)> 


(4tt) 


Wifi 


dXpdypdXpdyp 


rT  rR  rT  rR 


*4  ^ 4 

f n(rp-rp,rT) 


x CR(rp-rp)  PveTp)  p(eRp)  p*(eTp«)  p*(eRp.) 
X Cs  (rp-rp')  m(t+r-rp)  m*(t-r'p) 

X exp[-i  ap(t+r-rp)+  i tog  a'p(t+r-r'p)l 

^ n^P”rP,rT^  = ^^rP>rx)> 

CR(rp-rp)  = <R(rp)  R(rp) 

Cs(rp-rp)  = <(V(Z_T?)  * (cpT^PR^  1 lV'(z_T?  ) * Cepr^pR)  ] 


(59a) 

(59b) 

(59c) 


X exp[-i  kg  T)vz  + i kQ  rj'  ^]>  (59d) 

where  homogeneous  correlations  are  assumed.  In  this  expression  the  correlation  functions 
Cf|.  CR  and  Cg  are  for  the  shading  function,  the  reflection  coefficient  and  the  geometric 
scattering  coefficient,  respectively,  and  the  expectation  operator  separates  due  to  the  inde- 
pendences of  II,  R and  tj.  Since  the  correlations  vanish  for  I rp-rp  I » 0 this  integral 
can  be  reduced  to  a single  surface  integral  and  evaluated  for  the  cases  of  scattering  from 
the  interior  and  exterior  regions  separately  as  in  the  preceding  section. 

The  geometric  scattering  coefficient  correlation  function  can  be  approximated, 
where  vz  is  much  more  slowly  varying  than  tj. 


Cgfrp-rp)  - <[  (epp'H?)>R)  • (k  - (-i  kg  ^Z)“V)  exp(-i  kg  t)vz)  ] 

X [(epT+$pR)  ■ (k  - (-i  kg  I'p-y')  exp(i  kg  v'z  r?' ) ] > 

= { vzvz  - vz«  k0  "P"1  l ^FT^PR)  • V'  1 
- Vz(-i  kg  u?rl  [ (c'pj-^pr)  • V ] 

+ (-i  kg  vzyx  (i  kg  i^r1  [ (^PT-^PR)  •Vlt%T+'e'pR)  • V 1 
X Q2(kg  Vz,  kg  v',  rp-rp)  (60) 
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where  the  function 

CtykO  vz,  kQ  v'z,  rp-rp)  = <exp[-i  kg  r\vz  + i kg  r{v'z]  > 

(61) 

is  the  joint  characteristic  function  of  the  bottom  elevations. 

Assuming  a jointly  normal  distribution  for  the  surface  the  joint  characteristic  func- 
tion becomes 

Q2(k0  »z,  k0  vz,  rp-rp)  = exp[--=jK  (k0  hi>z)2  + (kQ  h^)2 

-2  (kQ  hpz)(kQ  hu')  Cv(  Irp-rp  1))  ] 

(62a) 

CT ?(  Irp-rp  1)  = <T?(rp)  rj(rp) > /h2 

(62b) 

n 

-a 

/— s 

o 

'w' 

II 

(62c) 

where  an  homogeneous  isotropic  bottom  is  assumed. 

The  above  expression  can  be  approximated  by  a first  order  expansion  of  the  sur- 
face correlation  function 

Cjj(  Irp-rp  1)  - 1 1 - Irp-rp  I/a2] , Irp-r^  1 < a 

s 0,  Irp-rp  1 > a , 

•(63) 

i.e.,  the  surface  correlation  vanishes  within  the  distance  a where  a is  the  surface  correlation 
distance.  Since  the  correlation  distance  is  short  compared  to  the  ranges  from  which  energy 
is  scattered,  the  angles  to  the  two  points  rp  and  rj>  can  be  equated. 

Collecting  results  then 

Cs(Vn>)  = | \ + 2vz(i  k0  vzY^  [ (ePT+epR)  *vl 

■ (ko  vz>~2  f (epT+ePR)  ' vi 2 | 

X exp  [-(kg  hvz)2  Irp-rp  l2/a2] 

(64) 

where  V'  - -V.  This  expression  can  be  readily  evaluated  with  the  substitution 

Irp-rp  I2  = (xp  - Xp)2  + (yp-yj>)2 

(65) 
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and  the  result  is 


Cs(rp-rp)  = | [v\  + i 4(h/a)2  v2  kQ[  (xp-xp)  t>x  + (yp-yp)  vy) 

- 4(h/a)4  v\  kj  [ (xp-xp)^x  + (yp-yp)^] 2 
+ 2(h/a)2  (^+„2)) 

X exp[  (kQ  h vja)2  ( (xp-xp)2  + (yP-yP)2)  ] (66a) 

vx  = (%T+/&PR)*t  (66b) 

vy  = (epT + %>r)  4 (66c) 

which  is  approximately  correct  for  all  xp,  yp,  xp,  yp. 

Collecting  results  the  correlation  function  of  the  reflected  field  is  given  by 


k02  rrrr  dxpdVpdxpdyp  „ , 

C(t  + T'  f)  ■ ~~i  JJJI  ~ C7r(rp-rp,rT)  C R(rp-rp) 

(47r)“  jjjj  rT  rn  rT  rn 


dxpdypdxpdyp 
rT  rR  rT  rR 


X I pC&j-p)  I 2 I p(^Rp)  I 2 m(t  + r-rp)  m *(t-r|j) 

X exp  (i  Wq  ap  (t  p-r  j>-r)  ) 
x j*£  + i 4(h/a)“i>“  kQ  [(xp-xp)i>x  + (yp-yp)^yl 
- 4(h/a)4  v 2 kg[  (xp-xp)t^x  + (yp-yp)^y  1 2 
+ 2(h/a)2  (»>2  + y2)} 

X exp[-(kg  hr»z/a)2  ( (xp-xp)2  + (yP-yP)2) ) * 

(67) 

The  integration  over  the  primed  coordinates  can  be  performed  after  expanding  rp  about 

Tp. 

From  equation  36e  then 

r p - rp  = tp  - tp  (68a) 
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and  analogous  to  the  expansion  in  equation  45  is  obtained 


( x'p  yp  + \ 

/ ^pt+®pr)  *^dx+  f (®pt+^pr)  ■'  dy  /co 
xp  Jy?  / 

(68b) 

where  Gpy  and  £pR  are  the  unit  vectors  at  the  scattering  point  rp  that  are  directed  along 
the  rays  to  the  transmitter  and  receiver,  respectively. 

With  reference  to  figure  3 the  expansions  for  the  direction  cosines  of  the  rays 
about  the  point  rp  are 


+ ^pr)  * 1 

COS0yx  + COS0Rx 

(69a) 

(£pT  + ^PR^  * J = 

COS0jy  + COS0Ry 

(69b) 

n 

It 

Cu 

i 

COS0JZ  + COSOj^j 

(69c) 

(SpT  +^pr)  " • - cos(0Tx  + + cos(0Rx  + eRx) 


— ^®pt  "*’<®PR^  * * 


sin-  0Tx  + sin2  0Rx  ^ (x^-xp  (69d) 


(SpT  +^pR)  ' J = cos(0Ty  + eTy)  + cos(0Ry  + eRy) 

f^PT  +^PR^  " j 


sin2  0 


Ty 


rT 


+ sin2  0Ry 
rR 


(yp-yp) 


^PT  +'GPR)  ' k = cos(0-pz+eyz)  + cos^Rz+eRz) 

(epT  +^pR)  ’ k 

where  the  angular  subscripts  are  defined  as  in  equation  46. 


(690 
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J 


figure  3.  Ray  vectors  at  the  points  r^#  and  "rjj  . 
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Substituting  the  direction  cosine  expansions  into  the  time  delay  expressions  and 
integrating  then 

rp-Tp  = (x  p - Xp>x/c0 


(sin2  0Ty  sin20Rv\  , 


+ (y'p  - yp)^/^ 


fsin^^jy  sin^  0 


rR 


Ly  (y'p-yp)2/2c0  (70) 


which  is  the  required  expansion. 

The  correlation  function  becomes 


C(t+r ,t)  = ^ — ff  d*p  d*P  I p(eTP)  1 2 lp(eRP)l2 
(4n)2JJ  4 rR 


X m(t+r-rp)  m*(t-Tp)  exp(-i  o>q  orp  r) 
X fj  dexdey  C^e.rj)  ( R(e) 

X jv2  - i W3)2  vl  kolexl'x  + W1 
- 4(h/a)4  V 2 k5[exvx  + eyi>y] 2 
+ 2(h/a)2  (r2  + v2) j 
X exp  [ - (kQ  hi>z/a)2  (e2  + e2)  ] 


X exp 


(/  sin2  0Tv  sin2  ( 


'Rx 


sin2  0 


+ eyvy  + 


Ty 


sin2  0 


Ry  < 


rR 


(71a) 
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where 


^ ■+.  A A 

rp-rp  = e = exi  + 6yj  • 


(71b) 


If  it  is  assumed  that  the  correlation  functions  of  the  shading  function  and  reflection  co- 
efficient are  approximately  constant  over  the  correlation  distance  of  the  medium,  the  inner 
integrals  become  elementary  integrals  and  the  correlation  integral  becomes 


As  for  the  case  of  the  mean  field  calculation  this  integral  is  evaluated  for  the  cases  of 
scattering  from  the  interior  and  exterior  regions  separately. 

Specular  Scattering  Correlation 

For  the  case  of  near  field  scattering  the  energy  is  received  largely  from  the  spec- 
ular point,  for  which  j>x  = v„  = 0,  and  the  correlation  function  is  approximated 


where  all  functions  in  the  integrand  are  approximately  constant  over  the  interior  region 
and  the  quantity  A is  the  area  of  the  ensonified  surface.  A crude  approximation  of  the 
area  A.  which  is  roughly  the  area  of  the  interior  region  and  that  gives  the  correct  range 
dependence,  is  given  by 


X m(t+7-rp)  m*(t-rp)  exp  (-i  WQOpr) 


v -2,2.  -2,  2,  2.2  2 , 2.  , 
X v [i>Z  + v (v  +v  ) - 2(v  +v  ) ] 
z z z x y v x y ' 


X exp[-  (a/2  hvz)2  (\+Vy)  1 ' 


(72) 


X m(t+r-Tgp)  m*(t-r§p)  exp  (-i  cjq  o^p  r)  (73a) 


TSP  RSP 


A a 


2 

(rTSP  + rRSP^ 


(73b) 
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With  this  expression  equation  73a  becomes 


C(t+r  ,t) 


Cn(°,rTSp)  CR(0) 

— j 

16(rTSp  + rRSp) 


^j2  I P(^TSP> 1 2 lP^RSP)  |2 


X m(t+r-Tgp)  m*(t-Tgp)  exp  (-i  ojq«spt)  (73c) 

which  is.  easily  obtained  and  shows  the  functional  form  of  the  correlation. 

Reverberant  Scattering  Correlation 

For  the  case  of  scattering  from  the  exterior  region  we  represent  the  complex  en- 
velope product  by  the  integral 

m(t+r-Tp)  m*(t-Tp) 

- T 

= % / dt'  m(t'+T)  m*(t' ) pe(t-Tp-t' ) (74) 

where  the  function  pe  is  defined  in  equation  51b.  Substituting  into  the  correlation  func- 
tion in  equation  72, 


«t+T,0  = tL™  2 exp  (a/21 ,)’  x / 


dt1  m(t'  + r)  m*(t'  ) 


Xff  ^ P 2~^  'P^TP^*2  (P(^Rp)  |2Cn(0,rT) 
rT  rR 

x [1  +(  l^pj+'SpR  I -2  ( |^pj-F8pR  I 2-v\  ] 

X exp(-i  cjq  ftp  t ) pe(t-rp-t')  x exp[-{a/2  hi>z)^ 

I&PT+^pr  1^1 

(75a) 
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where 


vx  + vy*vz.  = ^PT+^PR  |2  ' (75b) 

By  the  same  argument  applied  in  the  mean  field  calculation  the  radial  coordinate  can  be 
integrated  to  obtain 


C(t+r,t) 


CR(0)  Cjj(t) 
16tt 


exp  (a/2h)2 


X 


X 


dt'  m(t'+r)  m*(t') 


d0y 


Pj  ^ 

2 o I P (%-p)  1“  I P C^Rp)  1“ 
vZ.  r„ 


[ 1 + ( lepj  + epR  I 2 - v~  )/v^ 


- 2(lepT+epRl2  - vh/t?]  exp(-i  co0  ap  r) 


X exp[-(a/2  hi>7)2  lep-p+epRl2]  (-dpj/dt*  ) 


(t,0y,t  ) 

(75c) 


where  the  integrand  is  evaluated  at  the  ellipse-point  (t,  Oj,  t1).  Also,  since  in  the  exterior 
region  the  shading  correlation  function  is  approximately  independent  of  Oj  and  t , it  is 
denoted  simply  as  a function  of  t. 

Two  useful  expressions  for  computing  the  integrand  are  the  ellipse  relation  and  the 
range  relation,  respectively. 


rT  + rR  = c0(t+tD-t’) 

(76a) 

Tj  + = R + (Dr-Dj)  k 

(76b) 

where  the  quantity  R is  the  horizontal  range  vector  from  the  transmitter  to  the  receiver 
and  the  quantities  DR  and  D-p  are  the  distances  to  the  bottom  from  the  transmitter  and 
receiver. 
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With  appropriate  transpositions  each  expression  is  squared  in  equation  76  and  com- 
bined to  obtain  the  range  expressions 


c0^t+tD-t'  )2  - [R2+(Dr-Dt)2] 

2|co(t+tD-t')2  - ' [R  ( Dj^-Dj)1c]  | 


(77a) 


cfi(t+tD-f)2  - (R-+(Dr-Dt)2] 

Tp  = ~T : ri  (77b 

R 2jc0(t+tD-t')  +5rp-  [R(DR-DT)ft]j 

The  first  relation  gives  the  transmitter  range  coordinate  to  the  ellipse  point.  Expressions 
are  obtained  for  the  receiver  range  and  polar  coordinates  to  the  ellipse  point  by  first  sub- 
stituting equation  77  into  the  ellipse  relation  in  equation  76a  to  obtain  the  polar  coordi- 
nate relation 


&RP  • [R+(DR-DT)k]  = 


%P  • [R+(Dr-Dt)£] 


2 cQ(t+tD-f  ) (R2-KDr-Dt)2] 
{c^tftD-f  )2  + [R2-KDR-DT)2]} 


2 cg(t+tj)-tr  ) ejp  • [R“KDj^-Dj)k] 
{c2(t+tD-f)2+[R2+(DR-DT)2)j 


(78a) 


and  then  substituting  into  equation  77b  to  obtain  the  range  relation 

_ {^(t+tp-t1  )2  -2  cQ(t+tD-t’ ) %-p  • [R+(Dr-D-] 

R U ^ 2|c0(t+tD-t* ) “ ^TP  [ R+(DR-D-p 

(78b) 

With  these  coordinate  relations  the  integrand  in  equation  75c  can  be  evaluated. 

The  derivative  of  the  transmitter  cylindrical  polar  coordinate  is  given  by  the  deriva- 
tive of  the  transmitter  range  coordinate 
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(79a) 


dpT  J_  dpf  J.  _d_  2 2 

PT  dt'  = 2 “dT  = 2 dt’  (rT+DT)  " rT 


and  the  range  derivative  expression  follows  by  direct  computation 


dry 

“dt7 


= 4c0[c0(t+tD-t' ) - gyp  • [R+(DR-Dy)1c]  ] 

\rTrR  dt  7 

X |cQ(t+tD-t'  )2  - (R2+(DR-Dy)2]]  _1 
X [c^(t+tD-t>-  2 CQ(t+tj)-t  ) gyp  • [R+(DR-Dy)k] 

+ [R2+(DR-Dy)2l  J"1  (79b) 

where  the  range  relations  in  equation  77f  and  equation  78b  were  used. 

The  scattering  correlation  function  contains  the  quantities 

lepT  + = 2[  1+^pj  • ^pr] 

= 2 + [r“  +rR-D2]/rTrR  (80a) 


vz  “ Dj/rj  + Dr/i-r  (80b) 

which  can  be  evaluated  using  the  preceding  range  relations.  However,  since  in  the  far  field 
these  functions  are  nearly  constant  in  0 j,  it  will  be  convenient  to  use  the  far  field  limits 

I^PT+  ^PR|2  - 4>  vz  - 0»  rT~rR  - °°- 

Substituting  into  the  correlation  function  and  considering,  for  simplicity,  the  case 
of  R » Dr  - Dy,  then 
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cocr(°)  cn(t)  M2  -> 

C(t+T,t)  = exp(a/2h)2 


X f dt'  m(t'+r)  m*(t' ) 

J0 


X J dejjlplejp)!2  I p(^RP)  1 2 exp(-i  cog  <*p  r) 


X [ 1 + (2/^z)4  - 2(2/^z)2]  [c0(t+tD-t'  ) - &pp  ' R) 

X fcg(t+tD-t'  >-R2l-*  [cjj(t+tD-f) 

-2  c0(t+tD-t'  ) ^p  • R + R2]"1 

X expMa/h^])^^  . (t 

Further  the  case  of  times  is  considered  much  greater  than  the  pulse  length  t+tp  » t 
and  the  integrals  in  the  correlation  function  uncouple 
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c0CR(0)Cn(t)  /a\  2 exp  (a/2h)2 
<(t+r,tl  47r  (h)  [cQ(t+tD)-R2] 


I 

f dt1  m(t  + r ) m*(t' ) 


.2n 


[ d0T^lp(eTP)l2  I p(6'Rp) 1 2 exp  (-i  co0 «p t) 


X [ 1 + (2 lvzf  - 2{2!vz)2]  [c0(t+tD)  - ejp  • R] 
X (cQ(t+tD)2  -2  c0(t+tD)  eTp  • R + R2}-1 


X exp  (-(a/h  jO2] 


(t,0-p) 


(82) 


These  approximations  reduce  the  correlation  function  to  a simple  form  from  which,  in  the 
following  section,  a Fourier  transformation  will  give  the  GIPS  of  the  reverberant  scattering. 


4.  TIME  EVOLVING  SPECTRUM 

GENERALIZED  INSTANTANEOUS  POWER  SPECTRUM  (GIPS) 

The  GIPS  is  the  time  varying  power  spectrum  of  a nonstationary  process.  It  is  re- 
lated to  the  nonstationary  correlation  function  of  the  process  by  the  Fourier  transform 

P(t,co)  = Re  J C(t+r,t)  exp  (i  cj  r)  dr  (83) 

— OO 


1 ^ 

and  this  result,  which  is  developed  elsewhere,  *■  is  the  extension  of  the  Wiener-Kintchine 
theorem  for  stationary  processes.  It  resolves  the  time  dependent  power  of  the  process 
into  time  varying  spectral  components  and  can  be  estimated  by  the  sliding  FFT  calcula- 
tion.12 


37 


Early  Arriving  Signal 


The  early  arriving  signal  is  a mixture  of  the  direct  path  signal  and  the  specular 
scattered  signal.  The  GIPS  of  this  signal  is  given  by  the  Fourier  transform  of  equation 
58  where  the  correlation  function  of  the  specular  scattering  is  given  in  equation  73c. 
The  result  of  this  calculation  is  approximated  by 


where  beam  pattern  and  path  length  differences  have  been  ignored  between  the  direct  path 
signal  and  bottom  scattered  signal  and  the  mean  specular  scattering  in  equation  50a  is 
assumed  to  be  negligible  due  to  the  exponential  function. 

This  expression  shows  that  the  GIPS  is  the  sum  of  two  functions.  The  first  is  the 
GIPS  of  the  direct  path  signal;  the  second  is  the  GIPS  of  the  direct  path  signal  which  has 
the  Lloyd  mirror  effect  modulation  impressed  upon  it.  The  amplitude  of  this  modulation 
depends  on  the  bottom  characteristics,  but  for  typical  conditions  can  be  expected  to  be 
of  the  order  of  unity,  i.e.,  the  Lloyd  mirror  modulation  produces  surface  variations  that 
are  approximately  50  percent  of  the  GIPS  values. 

Late  Arriving  Signal 

The  GIPS  for  the  late  arriving  signal  is  given  by  the  Fourier  transform  of  equation 
82.  Since  this  expression  is  a product  of  two  functions  of  the  transform  variable,  there 
is  obtained  from  the  the  convolution  theorem 


P(t,w)  = • IpCejR)!2  ip(SRT)l2  m(t)  exp  [i  (co-<o0aD)t]  M*(to-a>0  aD) 


X jl  +Cn(0,rTSP)  CR(0)(a/4h)2  exp[i  0J0t(aD-aSp) 


-i(co-tOQ  a<jp)  rSp] 


(84) 


(85a) 


T 

P(co)  = Fr  f dt 1 m(t'+r)  m*(t') 


(85b) 
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Q(t,oj) 


X (1  + (2/t>z)4  -2(2 lvz)2]  [c0(t+tD)  - ejp  • iC] 


X (cQ(t+tD)  -2  c0(t-tD)  ejp  * iC  + R2]-1 


X exp[-(a/h  i^)2] 

) (t,0j) 


(85c) 


where  the  asterisk  denotes  convolution  and  the  symbol  Fr  {•}  denotes  Fourier  transfor- 
mation wrt  the  “t”  variable.  The  first  Fourier  transform,  P(co),  is  the  signal  modulation 
power  spectaim,  while  the  second  is  readily  found  to  be 


Q(t,ou)  = 2n  I d Oj  ^Me-pp)!2  lp(eRp)l ^ 6(co-cjq  ap) 

X [ 1 + (2/r>z)4  - 2(2/*z)2]  [c0(t+tD)  - ejp  • R] 

X (cQ(t+tD)2  -2  c0(t+tD)  eTp  • R + R2]-1 

Xexp  [-(a/hy  )2]|  • (86) 

I (t,0T) 

Inspection  of  this  integral  shows  that  the  delta  function  contributes  only  for  those  angles 
satisfying  the  condition 


oj  - gjq  c*p  = 0 , (87a) 

i.e.,  the  movements  of  the  transmitter  and  receiver  platforms  generate  a mapping  between 
angle  and  frequency.  Denoting  the  direction  of  transmitter  and  receiver  travel  to  be  a-p 
and  respectively,  then 


c0 


COS(0y  - ttp) 


COS  (6  - 0!j^) 


CJ-OJq 


Wq 


(87b) 
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This  expression  is  valid  for  the  far  field  in  which  the  vertical  deflections  can  be  ignored 
and  the  transmitter  and  receiver  polar  angles  are  related  through  equation  78a. 

This  relation  shows  that  the  energy,  transmitted  into  the  angle  0-p  and  received 
from  the  angle  0R,  is  mapped  into  a limited  band  of  frequencies.  Further,  it  shows 
that  the  energy,  transmitted  and  received  via  directions  that  are  symmetrical  about  the 
platform  movements,  are  mapped  into  the  same  frequencies;  i.e.,  defining  the  symmetri- 
cal directions  to  be  0y  and  0R  then 


Consequently  equation  87b  will  have  two  sets  of  angular  roots  (0j,  0R)  for  each  frequency. 
The  times  of  arrival  corresponding  to  each  set  of  roots  must  satisfy  equation  78a  and, 
generally,  will  not  be  the  same.  However,  simple  geometric  arguments  can  be  constructed 
to  show  that  equation  87b  will  always  have  a double  root  system.  In  the  long  time  limit, 
t -►  °°,  equation  78a  yields  the  limit  0R  - 0j.  Equation  87b  shows  that  in  this  limit  the 
angles  corresponding  to  the  maximum  and  minimum  frequencies  are  separated  by  it 
radians.  The  axis  that  separates  these  angles  is  designated  as  the  ambiguity  axis.  The 
integral  in  equation  86  separates  into  two  terms  corresponding  to  the  contributions  from 
each  side  of  the  ambiguity  axis  at  each  frequency 


(88a) 


aR  ‘ °R  ~ 0R  _aR  ' 


(88b) 


Q(t,co)  = Q(1>(t,co)  + Q(2)(t,co) 

Q(k)  (t,co)  = 27r  Ip(ejp)!2  lp(eRp)l2 


(89a) 


X [1  + (2 lvz)A  - 2(2 lvz)2)  [cQ(t+tD)  - eyp  • R] 
x [cQ(t+tD)2  -2  c0(t+tD)  eTP  • R + R2]-' 

X exp[-(a/h  i»z)2]  . 


(89b) 


Collecting  terms  the  GIPS  in  equation  85a  becomes 


X P(ou)  * [Q(D  (t,cj)  + Q(2)  (t,to) 


(90) 
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which  shows  that  the  power  at  each  time  and  frequency  is  contributed  from  both  sides 
of  the  ambiguity  axis.  Each  contribution  is  amplitude  weighted  by  the  beam  pattern 
and  transmission  path  associated  with  the  respective  side  of  the  ambiguity  axis  and 
the  power  is  smoothed  in  frequency  via  a convolution  with  the  transmit  power  spec- 
trum. If  the  bandwidth  of  the  transmitted  signal  is  small  compared  to  the  band  of 
Doppler  shifted  frequencies  then  it  is  approximated 

P(c o)  = 2ir  8 (oj) 
and  the  GIPS  obtained 


P(t,co)  - cq  Gr(0)  G'j-j(t) 


exp  (a/2h)2 

[c^(t+tD)-R2l 


X [Q<1)(t,co)  + Q<2)(t,co)]/2  • 


(92) 


This  is  our  principal  representation  of  the  far  field  GIPS  which,  to  be  complete,  needs 
an  explicit  angle-frequency  mapping  relation. 

An  exact  angle-frequency  mapping  relation  for  the  far  field  can  be  obtained  from 
equation  78a  and  equation  87b  but  for  these  purposes  it  is  instructive  and  useful  to 
obtain  a longtime,  t - zeroth  order  approximation.  Substituting  this  approximtion 
of  equation  78a, 


cos(t>R-D  = cos(0T-D  , (93) 

which  is  exact  in  the  long  time  limit,  into  87b,  the  angle-frequency  mapping  obtained 
is 


cos(0T-f)  = — I 


A = 
B = 


B / cj-coq\  C 


AM  <o0 


A2- 


'CO-OOq' 

CUq 


2-i 


1/2 


[v“  + 2vtvr  cos(a-p-aR)  + v^l  ^2/c0 
[vT  cos(r-aj)  + v R cos(f-aR)  ] /cQ 


C - (v-p  sin(r-a-p)  + vR  sin(r-aR)  ] /cq 


(94a) 

(94b) 

(94c) 

(94d) 


where  T is  the  angle  to  the  range  vector  R and  where  the  constants  A,  B,  C are  strictly 
functions  of  the  track  parameters.  The  constant  A gives  the  maximum  range  of 
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fractional  Doppler  shift  in  frequency,  B gives  the  fractional  Doppler  shift  in  frequency 
due  to  the  relative  range  rate  and  C gives  the  fractional  Doppler  shift  due  to  the  rela- 
tive transversal  rate.  In  principle  these  constants  are  measurable,  since  A and  B can  be 
estimated  from  spectral  data  and  they  are  related  to  C through  the  relation 

A2  = B2  + C2  • (94e) 

Further,  the  angle-frequency  mapping  shows  explicitly  that  two  angles  are  mapped  into 
each  frequency  in  the  Doppler  band  coq(1  ± A);  i.e.,  the  spectral  density  at  any  fre- 
quency in  general  cannot  be  unambiguously  resolved  into  components  associated  with 
the  directions  of  transmission.  The  direction  of  the  ambiguity  axis  corresponds  to  the 
upper  and  lower  edges  of  the  Doppler  band  and  is  given  by 

0^  = r + cos-,(B/A)  • (940 

These  mapping  relations  can  be  derived  more  accurately  by  using  higher  order  approxi- 
mations of  equation  78a. 

SPECTRAL  VARIANCE 

The  variance  of  the  spectral  ensemble  provides  an  estimate  of  the  usefulness  of 
any  one  spectral  estimate  and  a bound  on  this  function  is  given  by  ^ 


A 2 P(t,w)  < Jj  drdr'  (e  £(t+r)  <£*(t)  - C(t+r,t)l 

X [ 0*(t+r ' ) 0(t)  - C*  (t+r',t)  ] j 

X exp[i  u)(t-t')]  (95) 

For  the  case  of  ocean  bottom  reverberation  the  signal  can  be  approximated  as  a Gaussian 
random  variable  and  using  the  Gaussian  moment  expression 
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E(x]x2x3x4)  = E[x jx2]  E[ x3x4] 

+E [ x i X3 ] E[x2x4] 

+E[xjx4]  E[x2x3] 

-2  E[?c ! ] E [ X21  E[ x3 ] E[ x4]  (96) 

the  variance  expression  can  be  approximated. 

For  the  real  signal  case  the  bound  readily  is  obtained 

A2  P(t,oo)  < < #2(t)> <1  M(a>-a>0 1 2 > + P2(t,cu) 

“2  < 0(t)>2  I < M(oj-gj0)  > I2 

< 2[<  02(t)  > <IM(co-oj0)I2  > 

-<0(t)>2  I < M(co-ouq)  > I2]  (97) 

where 

P^ffiou)  = El  Re  0(t)  exp(i  cut)  M*(cu-cuq)  I2 

<<02(t)>  < lM(cu-cu0)l 2 > • (98) 

Normalizing  wrt  the  GIPS  envelope 

Pg(t,cu)  = < 02( t)>  < lM(cu-cu0)l2  > 

the  coefficient  of  variation  obtained  is 
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A2  P(t,«)  ^ „ 

“ ^ 2 

P|(t  ,co) 


<0(t)  >2  I < M(o>-cj0)  > I2 


(100a) 


<02(O>  < IM(w-(o0)I2> 


which  has  the  approximate  bound 


(100b) 


for  both  the  early  received  signal  and  the  late  received  signal.  Loosely  interpreted,  this 
result  shows  that  the  most  any  one  spectral  estimate  is  expected  to  differ  from  the  GIPS 
is  approximately  40  percent  of  the  GIPS  value.  This  variation  is  comparable  to  the 
amplitude  of  the  Lloyd  mirror  effect  modulation  of  the  GIPS  in  equation  84,  which 
implies  that  the  Lloyd  mirror  effect  cannot  be  estimated  from  one  realization  of  the 
spectrum  for  the  case  of  a moderately  rough  bottom. 
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CONCLUSIONS 


The  GIPS  expressions  obtained  here  for  the  early  received  signal  and  the  later 
received  signal  are  our  results  of  principal  interest.  They  exhibit  relations  between 
transmitter  and  receiver  parameters  and  the  received  spectral  components  that  can  be 
exploited.  In  particular,  the  expressions  for  the  later  received  signal  in  equation  89  and 
equation  92  are  a generalization  of  Halley’s  original  work  and  show  how  the  transmitter 
and  receiver  beam  deployment  parameters  influence  the  later  received  signal.  These 
expressions  provide  insights  into  ocean  information  transfer  problems. 


45 


REFERENCES 


1.  O.  D.  Grace,  Acoustic  Ray  Wave  Synthesis.  NUC  TP  408,  June  1974. 

2.  O.D.  Grace,  Signal  Variation  in  a Random  Medium.  NUC  TP  474,  July  1975. 

3.  Carl  Eckart,  The  Scattering  of  Sound  from  the  Sea  Surface.  J.  Acoust.  Soc. 
Am.,  25,  566-70,  May  1953. 

4.  C.  W.  Horton,  Sr.  and  T.  G.  Muir,  Theoretical  Studies  on  the  Scattering  of 
Acoustic  Waves  from  a Rough  Surface.  J.  Acoust.  Soc.  Am.,  41_,  627-34, 
March  1967. 

5.  C.  W.  Horton,  Sr.,  S.  K.  Mitchell  and  G.  R.  Barnard,  Model  Studies  on  the 
Scattering  of  Acoustic  Waves  from  a Rough  Surface.  J.  Acoust.  Soc.  Am., 

41,  635-43,  March  1967. 

6.  E.  P.  Gulin,  Amplitude  and  Phase  Fluctuations  of  a Sound  Wave  Reflected 
from  a Statistically  Uneven  Surface.  Sov.  Phys.-Acoust.,  j5,  135-140,  October- 
December  1962. 

7.  Petr  Beckmann  and  Andre'  Spizzichino,  The  Scattering  of  Electromagnetic 
Waves  from  Rough  Surfaces.  Pergamon  Press,  New  York,  1963. 

8.  I . Fortuin,  Survey  of  Literature  on  Reflection  and  Scattering  of  Sound  Waves 
at  the  Sea  Surface.  J.  Acoust.  Soc.  Am.,  47,  1209-28,  May  1970. 

9.  Robert  Halley,  A Method  for  Passive  Determination  of  a Range  to  a Pinging 
Searchlight  Sonar.  NEL  TM  231,  January  1957. 

10.  Robert  Halley,  Notes  on  the  Analysis  of  Reverberation  Signals  from  a Remote 
Pinging  Searchlight  Sonar.  NEL  TM  261,  June  1957. 

1 1.  Robert  Halley,  An  Application  of  the  Analysis  of  Reverberation  from  a Re- 
mote Searchlight  Sonar.  NEL  TM  262,  September  1957. 

12.  O.  D.  Grace,  Time  Evolving  Spectra.  NUC  TP  509,  April  1976. 


46 


APPENDIX:  HUYGENS  PRINCIPLE 


It  is  shown  elsewhere*  that  the  acoustic  field  within  a volume  V containing  har- 
monic sources  and  bounded  by  the  surface  S (see  figure  A.l)  is  given  by  RWS  formu- 
lation 


<A(r)  = + 


f dv  0(f)  -rdh 
I \/ 


+ f ds  n • 10(1)  v- V0(i)  1 
J S 


X(f.r) 


(A.l) 


where  the  unit  vector  n is  directed  into  the  volume.  In  this  field  expression  the  first 
term  represents  the  direct  path  field  and  the  volume  integral  represents  the  field  due  to 
scattering  by  the  medium  inhomogeneities.  The  surface  integral  represents  the  field 
scattered  by  the  surface. 

In  the  following  a theorem  by  Baker  and  Copson**  is  extended  and  shows  that 
the  RWS  surface  integral  produces  a null  effect  within  V when  the  surface  field  is  re- 
placed by  the  direct  path  field;  i.c.,  Huygen’s  principle  is  demonstrated  showing  out- 
going waves  have  no  effect  at  a point  after  propagating  past  that  point. 

Consider  a rederivation  of  equation  A.l  with  the  modification  that  the  sources 
are  excluded  from  V by  a surface  S' (see  figure  A.l).  By  inspection 


<£(  r)  = 


47T 


dv  0(|)  7(i,r) 


+ 


+ 


n • 10(|)  V-  V 0<£)1 


X(£.r) 


(A. 2) 


* O.  D.  Grace,  Acoustic  Ray  Wave  Synthesis.  NUC  TP  408,  June  1974  and  O.  D.  Grace,  Signal  Variation  in  a 
Random  Medium.  NUC  TP  474,  July  1976 

**  B.  B.  Baker  and  E.  T.  Copson,  The  Mathematical  Theory  of  Huygen’s  Principle,  Oxford  University  Press, 
1939. 
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— 


Figure  A.l.  Acoustic  cavity. 


where  n is  directed  into  the  volume.  Letting  the  surface  S be  a spherical  shell  of 
radius  R tending  to  infinity  then 


IS  = /dw  R2[0(£)Vr  -Vr  tff>]  X(I,r)  (A.3) 

where  dw  is  the  solid  angle  subtended  by  an  element  of  the  sphere.  For  the  modified 
radiation  condition  * * 


r!!Lr21^)Vr  - VR  0(1)  ] X(V)  = 0 (A.4) 

this  term  vanishes  and  by  redefining  the  unit  normal  vector  to  be  directed  within  S'  the 
“exterior”  form  of  the  RWS  formulation  is  obtained,  as  opposed  to  the  “interior”  form 
given  by  equation  A.  1 


<t>(r) 


4tt 


dv  0(f)  7(|, r) 


J ( ds'  H . [0(f)  v-V  0(1)] 


X(S,r) 


(A.5) 


Huygen’s  principle  is  demonstrated  with  the  aid  of  the  interior  and  exterior  forms. 

Consider  the  cavity  in  figure  A.  2 in  which  the  field  point  and  sources  are  sep- 
arated by  the  imaginary  boundary  B which  divides  S into  the  surfaces  Sj  and  S2,  i.e., 
S = Sj  + St.  If  S and  B are  imaginary  constructions  within  the  field,  then  the  fields 
on  S and  B are  given  by  the  direct  path  field  and  the  field  at  r can  be  represented  by 
both  the  interior  and  exterior  forms,  respectively 


<*r)  = 4^ 


J dv  0q(|)  7(|, r) 


+ f dsn  • [0O(|)V- V 0O(I)1 
J S.+B 


X(£,r)  (A.6a) 
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Figure  A.2.  Constructed  surfaces. 
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dsn  • [0o(£)V-Wo(S>1 


X(*.r)  • (A. 6b) 


+ 


/ 


S^+B 


Equating  expressions  the  integrals  over  B are  found  to  cancel,  due  to  opposing  unit  vec- 
tors, and 


^dsn  • [0o(l)V-V  0O(I)1  X(lr)  = 0 • (A. 7) 

This  result  shows  that  the  surface  scattering  operator  of  the  RWS  formulation,  with  re- 
spect to  an  arbitrary  closed  surface  S,  produces  a null  effect  within  S when  acting  on  the 
incident  field,  i.e.,  that  Huygen’s  principle  is  satisfied. 
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